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PREDICTION  OF  OXIDATION  HALF-WAVE  POTENTIALS  BY  MNDO  CALCULATIONS 
INTRODUCTION 

We  have  been  interested  in  the  possibility  of  predicting  half-wave 
potentials  for  molecules  being  considered  for  electrochemical  studies  to 
determine  if  these  potentials  will  fall  within  the  electrochemical  window 
of  the  solvent.  Calculated  molecular  orbital  energies  are  directly  related 
to  molecular  ionization  energies  via  Koopmans*  theorem  (1),  and  ionization 
energies  should  play  a  major  if  not  predominant  role  in  deteTiming  half¬ 
wave  potentials  (2,3);  therefore,  it  appeared  logical  to  seek  a  correlation 
between  calculated  molecular  orbital  energies  and  the  measured  half-wave 
potentials.  Indeed,  such  correlations  have  previously  been  successful  with 
Huckel-type  calculations  (2,4).  We  wish  to  report  here  the  first  use  of 
the  newly-developed  MNDO  (Modified  Neglect  of  Diatomic  Overlap)  semiempiri- 
cal  technique  for  these  correlations. 

The  MNDO  method  (5)  is  a  semiempirical  molecular  orbital  technique  of 
a  higher  level  than  generally  has  been  available  in  the  past.  The  method 
is  parameterized  by  comparison  of  calculated  results  with  experimental 
results  for  a  set  of  reference  molecules.  It  has  been  shown  to  give 
reasonable  predictions  for  the  ionization  energies  (i .e. ,  the  negative  of 
the  calculated  highest  occupied  molecular  orbital  [HOMO]  energy)  of  a  wide 
variety  of  molecules  (mean  absolute  error  of  0.48eV)  even  though  parameter¬ 
ized  specifically  for  molecular  geometries  and  heats  of  formation  (6). 

This  accuracy  should  allow  reasonable  predictions  of  oxidation  half-wave 
potentials,  given  the  success  of  the  simpler  Huckel  method. 

Of  particular  interest  in  this  laboratory  is  the  electrochemical  be¬ 
havior  of  compounds  in  AlCl^N-fr-buty^pyridinium  chloride  electrolytes. 


n 

which  are  molten  at  or  near  room  temperature  (7,8).  These  melts  provide 
a  completely  anhydrous  and  aprotic  medium  in  which  organic  electrochemistry 
may  be  performed. 

RESULTS 

Table  I  presents  the  results  of  the  calculations  for  our  correlation 
set  of  molecules  compared  to  both  measured  ionization  energies  and  measured 
half-wave  potentials.  The  mean  difference  between  the  predicted  ionization 
energies  and  the  measured  is  0.60eV,  but  the  predicted  value  was  consistently 
too  large.  The  correlation  was  made  between  the  predicted  ionization  ener¬ 
gies  and  the  oxidation  half-wave  potentials  in  acetonitrile/O. 50M  NaClO^ 

vs.  a  Ag/0.1MAg+  electrode.  A  least  squares  fit,  shown  in  Figure  1,  gave 

the  following  equation  for  the  prediction  of  half-wave  pu^ential  from  the 

calculated  ionization  energy: 

E,  =  0.725  E.  -  4.95  [1] 

-I 

E^  is  the  half-wave  potential  in  volts  and  Ej  is  the  calculated  ionization 
potential  in  eV. 

Robinson  and  Osteryoung  have  shown  a  good  correlation  between  ioniza¬ 
tion  energy  and  half-wave  potential  for  several  aromatic  species  in 
AlCi^rN- (n-butyl) -pyridinium  chloride  melts  (9).  They  obtained  an  almost 

constant  difference  of  0.23V  between  the  half-wave  potentials  measured  in 
+ 3 

the  melt  vs.  Al/Al  and  those  measured  in  acetonitrile/O. 50M  NaClO^  vs. 
Ag/0.1MAg+.  (Potentials  in  the  chloroaluminate  melts  in  both  the  Robinson 
and  Osteryoung  study  in  the  present  work  are  reported  with  respect  to  an 
aluminum  reference  electrode  immersed  in  2:1  A1C1 ^:N- (N-butyl)pyridinium 
chloride.)  We  modified  Eq.  [1]  for  this  difference,  and  obtained  the 
following  equation  for  prediction  of  half-wave  potentials  in  the  melt: 
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TABLE  I 


CORRELATION  BETWEEN  IONIZATION  ENERGIES  AND  HALF-WAVE  POTENTIALS 


Compound 

Calculated 

Ionization 

Energy— 

(eV) 

Observed 

Ionization 

Energy 

(eV) 

Half-Wave  Potentials 
in  Acetonitrile 
Relative  .  to 
Ag/O.lMAg  -  (V) 

(1) 

Benzene 

9.39 

9.25^ 

2.00 

(2) 

Toluene 

9.28 

8.82- 

1.98 

(3) 

o-Xylene 

9.23 

8.45^- 

1.57 

(4) 

m-Xylene 

9.24 

8.50^ 

1.58 

(5) 

p-Xylene 

8.73 

8.37- 

1.56 

(6) 

Biphenyl 

9.12 

8.20^- 

1.48 

(7) 

Naphthalene 

8.58 

8.15^ 

1.34 

(8) 

Anthracene 

8.09 

7.15^- 

0.84 

^MNDO  method,  present  work 
^Reference  6 
^Reference  11 
^Reference  12 
^Reference  13 


f 

-Reference  14 


CALCULATED  IONIZATION  ENERGY  (eV) 


[2] 


=  0.725  Ej  -  4.72  [2] 

E  is  the  half-wave  potential  in  the  melt  in  volts  relative  to  A1/A1+^. 

-m 

Use  of  these  equations  resulted  in  predictions  for  half-wave  potentials 
of  the  aromatic  molecules  given  in  Table  II.  These  predictions  also  were 
compared  with  Robinson  and  Osteryoung's  measured  values  in  the  melt  (9)  for 
three  species;  here  the  mean  absolute  difference  was  0.08V. 

Since  the  upper  end  of  the  electrochemical  window  of  the  room  tempera¬ 
ture  melt  is  only  about  2V  vs.  A1/A1+3,  the  predictions  preclude  the 
possibility  of  anodic  oxidation  of  the  nitroaromatics  in  the  melt.  The  pre¬ 
diction  for  thianthrene  is  for  a  very  low  oxidation  potential,  which  might 
indicate  the  possibility  of  spontaneous  oxidation  if  a  suitable  oxidant 
(or  electron  acceptor)  is  present  in  the  system.  The  electronic  and  ESR 
spectra  of  thianthrene  dissolved  in  an  acidic  AlCl3:N-(n-butyl)pyridinium 
chloride  solvent  suggest  that  the  radical  cation  of  thianthrene  is  formed 
spontaneously  (10).  The  measured  half-wave  potential  for  thianthrene  was 
1.17V  in  the  melt  (10),  about  IV  greater  than  the  predicted  potential. 

This  discrepancy  was  the  result  of  the  different  nature  of  the  HOMO  in 
thianthrene  compared  with  the  other  molecules.  In  the  case  of  thianthrene, 
the  HOMO  was  a  lone  pair  orbital  on  the  sulfur  atoms,  while  for  the  cor¬ 
relation  set  molecules  the  HOMO  was  a  tt  orbital. 


CONCLUSION 

In  summary,  the  MNDO  method  has  been  shown  to  reliably  estimate  half¬ 
wave  potentials  for  a  variety  of  aromatic  molecules,  providing  that  the 
same  type  of  HOMO  is  involved.  The  predictions  indicate  that  some  compounds 
may  be  oxidized  spontaneously  while  others  will  fall  outside  of  the  elec¬ 
trochemical  window  of  the  room  temperature  melt.  The  method  will  be  used 
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TABLE  II 


PREDICTED  HALF-WAVE  POTENTIALS  FOR  SELECTED  AROMATIC  MOLECULES 


Compound 

Calculated 

Ionization 

Energy 

(eV) 

Predicted  Half- 
Wave  Potential 
in  Acetonitrile 
Relative  to 
Ag/O.lMAg  V} 

Predicted 
Half-Wave 
Potential  in 
Melt  Relative 
to  Al/Al  -(V) 

Measured 
Half-Wave 
Potential  in 
Melt  Relative 
to  Al/Al  (V) 

Biphenyl 

9.12 

1.66 

1.89 

1.73- 

Naptha lene 

8.58 

1.27 

1.50 

1.58- 

Anthracene 

8.09 

0.92 

1.15 

1.1<£ 

Nitrobenzene 

10.31 

2.52 

2.75 

- 

m-Dinitrobenzene 

11.20 

3.17 

3.40 

- 

2,4, 6-Trinitrobenzene 

12.07 

3.80 

4.03 

- 

2,4,6-Trinitrotoluene 

11.83 

3.61 

3.84 

- 

Thianthrene 

6.70 

-0.09 

0.14 

1.17— 

According  to  Eq.  I 1 ] 


According  to  Eq.  [2] . 


Referene  9. 


—  Reference  10. 


to  estimate  half-wave  potentials  for  a  group  of  yet  unsynthesized  aromatic 
molecules  being  considered  for  electrochemical  studies. 
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